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I. 



INTRODUCTION 



Classical noncoherent signal detection is generally under- 
stood to mean the detection of a sine wave with random phase 
or time delay in additive white Gaussian noise (WGN) . However, 
the extension of this problem to a general noncoherent problem 
involving the detection of ah arbitrary periodic signal with 
random time delay in additive WGN, has received little attention 
A case of particular practical interest involves the detection 
of a baseband square wave with random time delay in additive WGN 

In digital communication systems, the message to be trans- 
mitted is encoded into a sequence of binary digits. Typically, 
these digits represented by the logical states * 1 ’ and 'O' are 
transmitted by sending a suitably chosen set of pulses which 
are distorted during both transmission and reception. The 
effect of this distortion is that at the receiver, it is no 
longer possible to determine exactly which waveform was actually 
transmitted. We can model the transmission as well as the 
distortion introduced in the receiver itself, as random noise. 

In all the analysis carried out, it will be assumed that the 
random noise can be modeled as additive, white Gaussian. The 
problem then is one of deciding, on the basis of noisy obser- 
vations, whether the transmitted waveform corresponds to a 
logical * 1 ’ or a logical 'O'. 

In this thesis, the analysis and design of noncoherent 
receivers that optimally (and in some cases suboptimally ) 
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detect an arbitrarily shaped periodic waveform with random 
time delay, is carried out. The random time delay is assumed 
to have some known probability density function (p.d.f .) . 

Two types of p.d.f.'s are assumed for the random time delay A. 
Namely, a uniform p.d.f. and a (.parameter varying) non-uniform 
p.d.f. have been considered, and their effect on the receiver 
probability of error (P ) has been studied. 

The design of optimum (in minimum sense) receivers is 
based on the likelihood ratio and the assumption of low signal- 
to-noise ratio (SNR) operations. The design of suboptimum 
receivers is based on heuristic approaches that intuitively 
yield reasonably good performance. 

Chapter II of this thesis presents statistical communication 
theoretic principles that have been applied to the analysis, 
design, and performance evaluation of noncoherent receivers. 
Approximations are made that lead to a decision rule or equiva- 
lently a receiver structure. A justification and meaning of 
such approximations is presented in Appendix B. It must be 
pointed out, that at high SNR, most receivers designed strictly 
from heuristic considerations perform adequately. At low SNR 
however, the receiver design must be optimum so as to not 
further degrade marginal operating conditions. Chapter II 
addresses this issue. 

Chapter III is devoted to the analysis and design of sub- 
optimum receivers where simple (heuristic) techniques have 
been applied to the design process. Performance evaluations 
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have been carried out, and wherever possible, comparisons are 
made with optimum systems operating in similar environments. 

The results of Chapters II and III involving receiver P e 
are analyzed and interpreted via the use of tables and graphs 
in Chapter IV. 

The fifth and final chapter presents some general conclu- 
sions to be derived from the work carried out in preparation 
of this thesis, and some suggestions for future analysis in 
this general topic area are given. 
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II. NONCOHERENT RECEIVER ANALYSIS 



A. BASIC CONCEPT 

Detection of arbitrary periodic signals with random time 
delay is a problem of significant importance in radar and 
communications, which has not received a great deal of attention 
in the past. The most relevant documentation related to this 
problem appears in the radar literature, where the problem is 
formulated as noncoherent detection of a radio frequency sine 
wave burst in the presence of noise [Ref. 1]. 

In communication systems, signals are transmitted to a 
receiver via some medium connecting the transmitter to the 
receiver called the channel. Transmitted signals undergo 
distortion in the channel as well as in the receiver itself. In 
many cases this distortion is caused by physical processes 
which, because of their complexity must be modeled via statis- 
tical means, i.e., random variables and/or processes. As 
previously pointed out in the receiver itself, noise is 
unavoidably added to the signals to cause further distortion 
and uncertainty. It is often found in practice, however, 
that the uncertainty created by the noise in a receiver results 
in the introduction of signal uncertainties, modeled as random 
processes having known waveshapes but random parameters such 
as random amplitudes, frequencies and/or phases. The problem 
of detecting signals of random time delay has been partially 
considered in Ref. 2. 
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B. PROBABILITY DENSITY FUNCTION OF A 

The type of random signal parameter dealt with in this 
thesis is the phase or time delay of an arbitrary periodic 
waveform of period T. Generally, the period T of the signal 
is very much shorter than the time duration T q of the signal 
and therefore it is extremely difficult to predict the time 
delay at the receiver of such a signal. In such cases it is 
reasonable to model the time delay 'A', as random variable 
having some density function that reflects the degree of 
uncertainty that exists about A. 

A useful p.d.f. for the time delay [Ref. 3] is given by 

2tt A 

mcos— pj7~ 

- e TI 0 (my A i- I T/2 I f 2 ’ 1 * 

0 _< m _< 00 

where the independent parameter 'm' determines the degree of 
uncertainty about the time delay A, and I^fm) is a modified 
Bessel function of the first kind of order zero. 

For m = 0, we see that f ^ ( A ) = ^ for A | T/2 | , so that 
A is uniformly distributed over the period T. As m increases, 
the uncertainty of the random variable A decreases and the 
p.d.f. approaches the shape of a normal density. As m ->■ <», 
f (A) tends to an impulse so that the time delay uncertainty 
tends to zero and the modeled random waveform approaches that 
of a deterministic or completely known signal. The shape of 
the p.d.f. as a function of m is shown in Fig. 2.1. 

The derivation of the mean and variance of this density func- 
tion is presented in Appendix A. 
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PDF OF TIME DELAY 




Figure 2.1 PDF of Time Delay 
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C. COMPUTATION OF THE LIKELIHOOD RATIO 



In this section, principles of statistical communication 
theory are applied to the derivation of a decision rule which 
will lead to the design of a receiver that optimally detects (with 
minimum probability of error, P ) a periodic signal with random 
time delay, in the presence of additive white Gaussian noise (WGN) 
under the assumption of low signal to noise ratio conditions. 

We start the problem by considering two hypotheses H^ and 
Hq such that during the observation interval (0,T ), under 
hypothesis H^ , we assume that the periodic signal with unknown 
time delay is present and under the second hypothesis, i.e., 

Hq , there is no signal present. In both cases, the effect of 
the WGN must be considered. Thus symbolically. 



H. 



r(t) = v(t-A) + n(t) 



H. 



r ( t ) = n ( t ) 



0 < t < T 
— — o 



( 2 . 2 ) 



where 



r(t) = the signal at the front end of the receiver; 

v(t) = the T-periodic deterministic signal with 
A a random variable modeling the unknown 
time delay of the waveform; and 

n(t) = a sample function of a white Gaussian process 
with zero mean and two-sided power spectral 
density level N q /2 watts/Hz. 

In order to satisfy the minimum probability of error cri- 
terion, the optimum decision rule is obtained by comparing the 
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likelihood ratio (XR)_, with a threshold n and decide in favor of 
hypothesis if LR is greater than ri , or in favor of Hq if 
LR is less than n. 

Due to the random variable A, v(t-A) itself is a random 
process. If we condition on A, i.e., we assume some fixed 
value of A, then v(t-A) is completely known. However, the 
observed signal r(t), conditioned on A, because of the noise 
n(t), is a Gaussian random process. In other words condition- 
ing on A, under either hypothesis, r(t) is Gaussian so that 
the likelihood functions can be expressed in closed form and 
are given by [Ref. 4] 

T o 

L (r) = F exp{ - / r 2 (t)dt} (2.3) 

0 o 0 



and 

T o 

L.(r|A) = F exp {- / [r ( t) -v ( t-A ) ] 2 dt } (2.4) 

o 0 

The constant F is of no consequence here as it will 
cancel once the LR, namely L^frl/L^r) is evaluated. Observe 
that 



00 

(r ) = / L 1 (r | A) f A ( A) dA 

— oo 



(2.5) 
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The likelihood ratio test is therefore 



L 

L 



1 

0 



(r) 

Try 



H, 




H, 



( 2 . 6 ) 



where 



and 



P{H 0 } C 10 - C Q0 

mqT c 01 - c 
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P{H i > 



i = 1,0 



(2.7) 



( 2 . 8 ) 



are the prior probabilities of occurrence of hypotheses H^, 
and 



C.. i , j = 0 , 1 (2.9) 

are the costs associated with making decisions about hypotheses 
H^. For a minimum probability of error receiver 

| 1 i ^ i 

C. . = (2.10) 

1 0 i = j 

Taking advantage of some simplifications in writing the LR, it 
is simple to show that the LR becomes 
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oo 



r (t) v (t-A) dt} 



L 1 (r) 

L 0 ( I } 



/ exp{^~ j 



o 0 



exp{ - A J ° V 2 (t-A)dt)f (A)dA 
o 0 



H ] 

> 

< 

H, 



( 2 . 11 ) 



Since v(t) is T-per iodic, it can be expressed in terms of an 
exponential Fourier series, namely 



v (t) 





( 2 . 12 ) 



where 



V, 



1 , T/2 

— / v(t) e dt for all k 

1 -T/2 



(2.13) 



Using the Fourier series expansion, we have 



T o p oo co _jll(k + £)x T o j^(k+£)t 

/ v (t-A)dt = l l V V e / e dt 

0 k=-°° Z=-co 0 

(2.14) 



However note that 



T o j^(k+£)t 
/ e dt 

0 



( T if l = -k 

) ° 

^ 0 otherwise 



(2.15) 
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provided that T q = nT, where n is an integer. This result is 
also valid even if n is not an integer, provided that n >> 1. 
Thus applying this result to Eq. (2.14) results in 



/ ° v 2 (t-A)dt = T l V k V* = T I | V | 2 = E (2.16) 

0 k=-°° k=-°° 



Observe that E is the energy of the signal v(t) for 0 <_ t <_ T . 
Furthermore , 



/ r(t)v(t-A)dt 



.2tt, , T 

00 , o 



j^-kt 



l V k e 1 / r ( t) e dt 



k=-oo 



l W 

k=-oo 



• 2tt, . 



(2.17) 



where 



0 



T ^TTk. 

/ r(t)e dt 



(2.18) 



So the test of Equation (2.6) becomes 



.2TTk. H ] 

oo — J X 

/ exp l V k R k e T }f A (A)dA > H exp(E/N A > 



o k=-°° 



'A 



H, 



0 



(2.19) 
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D. TWO TERM APPROXIMATION 

In order to proceed with the analysis so as to derive a 
signal processing algorithm from Equation (2.19), we will 
assume that 



2 

rp- y V, R, e 
N , L k k 
o k=-°° 



. 2-rrk 

-D-sr- A 



<< 



( 2 . 20 ) 



In Appendix B, we will show under what conditions this 

assumption is valid. It turns out that the assumption of 

Equation (2.20) is essentially equivalent to an assumption of 

low SNR, i.e., E/N << 1. 

o 

Since 



e x - 1 + x if x<<1 



( 2 . 21 ) 



provided Equation (2.20) is valid, it is possible to approxi- 
mate the exponential term of Equation (2.19) by two terms only 
so that the test can be written as 



/ (1 



— I 

N , t- 
o k=-» 



\ R k e 



. 2-nrk 
-D-m-A 



) f A ( A ) dA 



H. 

> 

< 

H, 



exp{E/N Q } 



( 2 . 22 ) 



or equivalently 




W 



/ 



— oo 




(X)dX 



H- 



> 

< 




N 

-y-tn exp{E/N o )-l] 



A 

= n* 



(2.23) 
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The test of Equation (2.23) is valid provided Equation 
(2.20) holds. 

Observe that in order to limit the approximation error to 
10% in Equation (2.21), |x| <_ 0.39. 

In the remaining analysis we assume equal prior probabili- 
ties, namely P{H Q } = j = P{H^} and cost assignment = 1 - 6^ , 

so the threshold of Equation (2.23) becomes 



In order to further analyze the test of Equation (2.23), we 
consider the following two cases: 




1 ] 



(2.24) 



1. Uniform p.d.f. on A 



2. Nonuniform p.d.f. on A 



1 . Uniform p.d.f. 



For this case, the p.d.f. of A is 




1 

T 



A < I T/2 | 



(2.25) 



and the test becomes 




(2.26) 
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However, since 




-T/2 




j T if k = 0 

^ 0 otherwise 



(2.27) 



the test of Equation (2.26) reduces to 

*1 

V 0 R 0 < n * 

H n 



(2.28) 



where is just the d.c. component of v(t) . Note that for 
k = 0, Equation (2.18) becomes 



T 

o 

R 0 =I r ( t) dt (2.29) 



so that the test can be expressed as 



n*/V Q (2.30) 



and the receiver structure takes on the simple form shown in 
Fig. 2.2. 




Figure 2.2 Nonchoherent Receiver (Uniform p.d.f., 
2 Term Approximation) 



H, 



/ r (t) dt 



> 

< 

H, 
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It is clear from Equation (2.30) that this simple receiver 
can only operate provided Vq ^ 0.. If V Q is indeed zero, some 
other approach must be used for the purpose of deciding on 
which hypothesis is true upon reception of r(t). In the sequel 
we discuss other detection approaches which will work even if 

V 0 = °- 

The performance of this receiver is evaluated as follows. 

Let 



£ 




r(t) dt 



(2.31) 



Conditioning on A, we see that £ is a conditional Gaussian 
random variable, where 



E{£ | H q } 




(2.32) 



E{ £ | H 1 ,X} 




v ( t-A ) dt 






- v o Z v k e 

k=-°° 



T . 2 irk. 
o D— t 

J e dt 



= V^T 
0 o 




(2.33) 



Furthermore 
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VarU| H 1 ,A} 



= E { [ U | H x , X ) - E{Jl|H lf A}] 2 ) 



T T 

o o 9 

= E{ [V / [v ( t-A ) + n ( t) ] dt - V / v(t-A)dtp} 

u 0 u 0 



T 

o o N 

= V / / E{n ( t) n (t) }dt dT = -f- V. T 
o ' 20o 



= Var { Jt | H Q } = of (2.34) 

Equation (2.33) shows that the conditional mean m^ is indepen- 
dent of A and Equation (2.34) also shows that the conditional 
variance of Z is independent of A. 

Hence 



-U-m ) 2 

f_(£|H lf A) = f ( £ | H-. ) = — exp{ (2.35) 

2o Z 

Also, from Equations (2.32) and (2.34) 

1 Z 2 

f T U|H n ) = — - — exp{ (2.36) 

0 2 °l 

The p.d.f.'s of Equations (2.35) and (2.36) are plotted in 
Fig. 2.3. 
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Fig. 2.3 Probability Density Function of £ 



From Equations (2.35) and (2.36) and Fig. 2.3, it is clear 
that the two p.d.f.'s are completely identical with different 
means. In order to have minimum probability of error, the 
threshold rf should be the point where the two p.d.f.'s inter- 
sect. That is, 



n 




v n T 

0 o 



(2.37) 



This threshold differs from that given by Equation (2.24) . 
We must however keep in mind that the results were obtained by 
approximations on the LR test of Equation (2.19). 

Receiver probability of error , is thus given by 

P = P (H n } P + P{H.)P M (2.3 

e O F 1 M 



where 
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P„ = Probability of false alarm; 

r 

P^ = Probability of miss; and 

P(H^} , i = 0,1 are the prior probabilities. 

Since we have assumed equal prior probabilities, P^ becomes 



n 2 2 n + , n , 2 2 

1 00 1 l /2a . i i (Jl-nu) /2 a. 

y / — - — e ^dii + y / — e 1 ^dS, 



H* /2tt a. 



/Jia, 



1 

2 



r 1 -x /2 , i r 

J e 7 dx + y J 






n*/a £ /2? 



i -y / 2 , 

e 1 ' dy 

/2ir 



i n* o*- m i 

y[erfc* (— ) + erf* ( — — =■) ] 



(2.39) 



From Equations (2.24), (2.33) and (2.34), we get 



P 

e 



y{erfc* [- 



E/N E/N 

° -1, , ^ r e ° -1 

— - — ] + erf* [ 



' 2v n T 
0 o 

N 



2 V n T 
0 Q 

N 



0 O 
N 



]} 



(2.40) 



Observe that if v(t) has a very strong d.c. component in 

2 

comparison to its harmonics, then VqT q /N o is approximately 
equal to the SNR defined as E/N . Then Equation (2.40) becomes 
a function of SNR only. Otherwise, it is a function of both 
the SNR and the V^T /N ratio. 
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2 . 



Nonuniform p.d.f. 



For this case, the p.d.f. of A as given by Equation 



(2.1) is 



f 



A 



(A) 



mcos 

e 



2tt A 

T 



TI 0 (m) 



A £ | T/2 | , 0 <_ m <_ 00 



It is known that [Ref. 5] 



m cos 0 
e 



I 0 (m) 



+ 2 J I (m) cos p0 
p=l P 



00 

1 I (m)cos p0 
p=-oo P 



(2.41) 



since 



I 

P 



( •) 



I 

"P 



(•) 



Therefore, Equation (2.1) becomes 



f A m 



y I (m) cos p 

u D 

p = -oo ^ 

TI Q (m) 



2 7T A 



(2.42) 



From Equation (2.23) we have the test 



. 27Tk, 

^ 00 “J-m— * 

I V k R k / t a) e dA 

K= — 00 — 00 



H ] 

> 

< 

H, 



n* 
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or equivalently 



where 




H. 




H, 



(2.43) 



F (^) 
A 1 T ' 



GO 




f A (A) e dX 



(2.44) 



Substituting Equation (2.42) into (2.44), we get 



F (^S) 
T 1 



TI n (m) ^ p 

0 p= — oo r 



co T/2 

l * / cos P 



-i— X 

2ttX 3 T A 



dX 



-T/2 



” T/2 j^(p-k)X 

£ I (m) [ / e 1 dX 



2TIq (m) P -t/2 



T/2 -j^L(p+k)X 
+ / e dX] 

-T/2 



(2.45) 



Note that both integrals in Equation (2.45) are zero unless 
p = ±k. For p = k, first integral in brackets is T, and the 
second integral becomes zero. For p = -k, the first integral 
is zero but the second integral is T. 
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Therefore 



,27rk, _ I k 

A ' T I 0 (m) ' 



(2.46) 



since 



W 



) = 



V). 



Hence the test of Equation (2.43) becomes 



__1 

Iq (m) 



l V k R k I k (m) 




H, 



(2.47) 



The test can be simplified further by letting 



I 0 (m) k=- V k R k I k (m) 



-| U w VJ 

rw {v o I o tm) / r(t) + I Vk (, i / r(t) 

1 0 ^m; 0 k=1 K 0 



. 2iTk. 

D— t 
e 1 dt 



T .2iTk, 

°° , ° l-m -1 

+ l V,I,(m) / r ( t ) e dt} 

-k=l K K 0 



_1 

I 0 (m) 



( v 0 I o (m) 



T 

o 

f r(t)dt 

0 



oo 

+ 2 l |v |I (m) 

k=l K K 



T o 

/ r(t)cos(^t + a k )dt} 



(2.48) 
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where 



V 



±k 




(2.49) 



Finally the test becomes 



T o co T o H 1 

V q / r(t)dt + j— jcrr l |v |l (m) / r ( t ) cos (^t + a ) dt J n* 

0 o V m) k=i k k o T * „ 

H o 

(2.50) 



and the receiver may be implemented as shown in Fig. 2.4. 




Figure 2.4 Noncoherent Receiver (Nonuniform p.d.f., 
2 Term Approximation) 
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Notice that this receiver can operate on signals both 
with and without d.c. component. When a signal has zero d.c. 
term, the receiver of Fig. 2.4 can be simplified by eliminating 
the uppermost branch. 

The performance of this receiver is obtained by using 
the fact that conditioned on X, under either hypothesis, g 
is a Gaussian random variable with 

E{g | H q } = 0 (2.51) 



and 



E{g| A} 



E{ 



1 

IqW 




W * 1 



T . 2tt k , 

, O D-m-t . 

/ v(t-A)e dt} 

0 



1 0 (m) 



l l 

k=-oo £ = -c 



v kVk (m)e 






T o j4fu+k)t 
/ e dt 

0 



(2.52) 



The integral of Equation (2.52) has been analyzed in Equation 
(2.15). Therefore 



E{g|H 1 ,X} 



T 

o 



I Q (m) 





(m) 




= h ( A ) 



(2.53) 



Furthermore , 
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Var{ g | H x , A } 



1 



oo 



„ T 27Tk. 

1 “ o D-m-t 

= E{ (y Try l V 1 (m) / [v ( t-A ) +n ( t) ] e 1 dt 

■ L 0 ' m ' k=-°° K K 0 



T . 2 ltk. 

O D-m-t 



k= L v k i k (m) j v(t - x)e T dt > 2 > 



T . 2irk t 

El ‘i^TST JLv k <-> J ° n(t,e T dt)2> 



OO oo 



- — I l V k V s, I k (m) I £ (m) // E{n (t) n (t ) } 



lQ(m) k=-°° £=- c 



.2irk, .2 tt£ 

D-m-t D Trn T 

e e 1 dt dx 



N oo oo T o j^(k+£)t 

l l V,V I (m)I (m) / e dt 



Iq (in) 2 k=-co £=-c 



0 



T N 

o o 



l |vJ 2 iS(»,) 



2 Iq ( m) k=-» k k 



= Var {g | H^} = Var{g|H Q } = 



A 2 



(2.54) 



With 



p = P {Hq } and 1-p = P{H^}, we have an expression for the 



conditional P g given by 
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oo 



P e (A) = p / f Q ( g | H Q ) dg + (1-p) / f ( ,(g|H 1/ X)dg 



= P f 



2 2 

1 ~9 /2a 



H * /2 tt a 



dg + (1-p) / 



n * 1 -U-h(\)) 2 /2ol 

/ 2tto 



dg 



n* n*“h(A) 

p erf c *(— ) + (1-p) erf * ( - ) 

9 ■ 9 



(2.55) 



where n*> h(A) and o are given by Equations (2.24), (2.53) 

9 

and (2.54), respectively. 

The actual receiver is obtained by integrating 
P (A) over the p.d.f. of A, resulting in 



P o = p erfc*(— ) + (1-p) / f A (A) erf * ( 



n*-h ( A) 



-) dA (2.56) 



Observe that 



n* 

a 



N E/N 0 

4“(e -1) 



1 1/2 



N T 
o o 



21 q (m) k=-«> k 



l |V v | 2 I^(m) 



E/N 



- 1 



W- . I l\l 2 l k ( ”>/ l 2 o< m ) ( l l v < 

o k=-°° / x.=-°° 



1/2 



(2.57) 



Also 
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h(A) 



a 



g 





(2.58) 



E. THREE TERM APPROXIMATION 

Whenever Equation (2.20) is not satisfied to the degree 
that Equation (2.21) is valid, it is possible to proceed with 
the analysis beyond Equation (2.19), if in Equation (2.19) we 
expand the exponential term and approximate it with three 
terms. Since 



in order to have the approximation error not exceed 10%, we 
must have x _< 0.79. Using the three term approximation of 
Equation (2.59) and applying it to Equation (2.19) , the test 
then becomes 



exp{x} 



2 

1 + x + h 



if x << 1 



(2.59) 



— CO 




(2.60) 
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or equivalently 




f A ( A ) dA > y 



(2.61) 



where 



E/N 



o 



Y 



- 1 ) 



(2.62) 



Here again, two cases must be considered: 

1. Uniform p.d.f. on A 

2. Nonuniform p.d.f. on A 
1 . Uniform p.d.f. 

For this case the first term on the left hand side (L.H.S.) 
of the inequality in Equation (2.61) has already been analyzed 
in Equations (2.26) through (2.28). This first term is just 
VqR q . The second term of Equation (2.61), i.e., 




becomes (with the aid of Equation (2.15)) 



A _J_ 

2T * 2 
N 



o 



k=-oo £= 



l 




i i\i 2 'V 2 

; — CO 



(2.63) 
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The test therefore becomes 



V 0 R 0 

N 



oo 

4 i 

N k=- c 



V, 



H ] 

> 

< 

H, 



1 = y- 

2 Y 



12.64) 



or equivalently (using Equations (2.13) and (2.18)) 



T 

V Q o 



T 

V„ o 



N / r (t ) dt + / r ( t) dt] 



o 0 



o 0 



H, 

oo I V I T 2 

+ 2 l [-*p- f ° r ( t ) cos ^^tdt] * y* 
k=1 ° 0 H q 



(2 . 65) 



This test leads to the receiver structure shown in Fig. 2.5. 

Note that this receiver utilizes both the d.c. and non- 
d.c. components of the periodic signal. However, in practice 
there are many cases where the d.c. component of the signal is zero. 
For those cases, the test reduces to 

H i 

4 i i v vi 2 i R ki 2 < y * <2 - 66) 

N C Hq 

and the corresponding receiver can be implemented by simply 
eliminating the upper most branch of the receiver in Fig. 2.5. 

2 . Nonuniform p.d.f. 

Our starting point here is Equation (2.11), repeated 
here for convenience. 
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Fig. 2.5 Noncoherent Receiver (Uniform p.d.f., 
3 Term Approximation) 
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00 9 ° 1 ° 9 

/ exp{— / r (t)v(t-A)dt}exp{- — / v (t-A ) dt } f ^ (A ) 



dA 



o 0 



o 0 



H ] 

> 

< 

H, 



where n is given by Equation (2.7). Note that the second 
exponential of the inequality given by Equation (2.11) has 
already been analyzed in Equations (2.14) through (2.16). 
Furthermore, if we expand the first exponential in the above 
inequality and approximate it with three terms as given by 
Equation (2.59), then the test of Equation (2.11) becomes 



oo 



/ 



— 00 





r(t)v(t-A)dt f^(A)dA 



+ 



1 

2 





r (t) v ( t-A ) dt] 2 f^ ( A ) dA 




(2.67) 



where y is given by Equation (2.62). 

Interchanging the order of integration in Equation 
(2.67), we obtain 



jj— / r (t) / v(t-A) f A (A)dAdt 



o 0 



A 



H. 



+ — 5 - // r (t) r (x) 

N 0 
o 



/ v (t-A) v (x-A) f^ ( A) dA dt dx < y 



( 2 . 68 ; 
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Note that v(t) is a deterministic signal but as mentioned 
earlier, due to the r.v. A, 

A 

y ( t) = v(t-A) (2.69) 

itself is a random process with expected value 



E{y(t)} = / v (t-A) f . (A ) dA = m (t) (2.70) 

— oo 

and autocorrelation function 

OO 

E{y (t ) y ( t ) } = / v(t-A)v(x-A)f A (A)dA = R y (t,x) (2.71) 

— OO 

Substitution of Equations (2.70) and (2.71) 
in Equation (2.68) leads to the following test that a receiver 
has to perform 

T T T H 1 

2 ° ? o o 

— / r(t)m (t)dt + —j / r(t)[ / r (x ) R ( t , x ) ] dt y (2.72) 

o 0 Y N O 0 Y „ 

° H 0 

and the implementation of this receiver is shown in Fig. 2.6. 
Here the autocorrelation function R (t,x) can be considered 
as the impulse response of a filter that has to be designed 
for detection of the signal. 

Furthermore, for the uniform p.d.f. on A, the expected 
value of y(t) given by Equation (2.70) becomes 
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Fig. 2.6 Noncoherent Receiver (Nonuniform p.d.f., 
3 Term Approximation) 



T/2 

itu.(t) = =• / v (t-A) dX 

x i _ T/2 



T/2 -jlliix 

i T V. e T / e T 
T k=-°° k -T/2 



dA 



= V, 



(2.73) 



which is a constant, and the autocorrelation function of y(t) 
given by Equation (2.71) turns out to be 
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1 

T 




-T/2 



T/2 “ j~7fr (k+£ ) A 

/ e dA 




00 



I l v , 




(2.74) 



which is a function of (t-x) only. This implies that for the 
uniform p.d.f. on A, the impulse response of the filter turns 
out to be time invariant. That is, it depends only on the 
time difference (t-x). Furthermore, m^(t) being a constant 
can be eliminated from the upper branch of the receiver shown 
in Fig. 2.6. 

F. ANALYSIS OF TWO SPECIFIC SIGNAL WAVESHAPES 

It was pointed out earlier that the receiver derived using 
a two term approximation on the exponential appearing in Equa- 
tion (2.19), is unable to discriminate between the signal 
present versus the signal absent case for signals having zero 
d.c. component when the random time delay obeys a uniform 
distribution. When the random time delay obeys a nonuniform 
p.d.f., a receiver was derived (see Fig. 2.4) that could 
detect the presence of signals having zero d.c. component. 

In this section, assuming a nonuniform p.d.f. on A, the 
receiver performance will be analyzed for two different signals 
which are found quite often in practice and that have a zero 
d.c. component. The two signals are 
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1 



Sine Wave 



2. Square Wave 
1 . Sine Wave 

For the first case let 



v ( t ). 



A cos 







(2.75) 



This periodic signal has only two discrete components. 
In other words, it can be represented in terms of its two 
Fourier coefficients, i.e., 



V 



k 




k = ±1 
otherwise 



Thus , 



(2.76) 



I K 



k=-oo 



<f > 2 ♦ <f , 2 



(2.77) 



and from Equations (2.57) and (2.58), we obtain 



n * 
a 



E/N 

, o . . 
(e -1) 



If (m) 1/2 

T - v , (2E/N ) L/ 

I ^ (m) ' o 



(2.78) 
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cos 



2tt A 

T ' 



(2.79) 



n *- h (.A ) 



, E / N o ,> , 2E I l (m > 
(e 

I 1 ^ 1/2 
ltw (2E/N o> V 



Thus the expression for becomes 



= p erfc* 



(ne 



E/N 



- 1 ) 



Il(m) 3/2 

rTHT (2E / N o» 7 

0 



(1-p) T/2 

+ TI (m) / erf * 

U -T/2 



(ne 



E/N 



o N 2E I 1 (m) 2 tt A 

- 1 *- FT - T^imT cos ~T~ 

I 1 1/2 

TJ5T (2E/ V 7 



2 tt 

mcos— A 
e 1 dA 

(2.80) 



Computation of Equation (2.80) can be carried out as 
a function of the SNR E/N q , p and m. This has been carried out, 
and the results are presented graphically in Chapter IV. For 
the sake of computational simplicity, equal prior probabilities 
were assumed, i.e.. 



P = 



1-p = 



1 

2 



resulting in n = 1 • 

2 . Square Wave 

For the second case, let v(t) be a periodic square 
wave as shown in Fig. 2.7. 
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Fig. 2.7 Periodic Square Wave 



The Fourier coefficients of v(t) are 



V, 



1 T/2 

mi j Ae dt 

0 



T/2 



T 

/ Ae dt] 



rk _ 

in* Ac -^ k r sin 2 i 

2 Ae 1 irk J 



Observe that V q and all the even coefficients are zero, i 



V 2k = 0 for all k 



and 



V 



±k 




2A 
kn ' 




,2A, 2 
kir 



/ 



k odd 



(2.81) 



S • ^ 

(2.82) 



(2.83) 
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Note that 



T/2 2 2 00 2 

f v ( t) dt = A T = T l |vj 
-T/2 k=-°o K 



which implies that 




Furthermore 





2l k (m) 



|v o | 2 l o <m) + 2 ^ l v kl 2 l k (m) 



so that 



oo 




I k (m) e 






I v o S 2 1 o (n>) 



+ 2 




2 l k (m) 



COS 



T A 



since 



|V k | = |v_ k I ; I k (m) = I_ k (m) 



(2.84) 



(2 . 85) 



( 2 . 86 ) 



(2.87) 
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Hence Equation (2.57) becomes 



n* 

a 



E/N 

(ne °-l) 

4E/N o k I 1 l V k| 2l S (m) 
A 2 Iq ( m) 



( 2 . 88 ) 



and Equation (2.58) becomes 



N 4 E/N oo 

n.-h(A) I k‘"» 




(2.89) 



Thus Equation (2.56) for the receiver probability of error 
becomes 



P = p erfc* 
e 



E/N 

/ o , . 

(ne -1) 



4E/N o l l V k! 2l k (m) 
k=l 

2 2 
Ig(m)A Z 



1/2 



N 



4E/N 



( ne "°-l) ( l |v k | 2 I k (m) oos^X) 



+(l-p) / erf* 
0 



I Q (m)A k=l 



4E 

N 



o k=l 



l l V k | 2l k (m) 



1/2 



7 2 

I 0 (m) A^ 



me os 



2tt A 



TI 0 (m) 



dX 



(2.90) 
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Observe that Equation (2.90) (as well as Equations 
(2.57) and (2.58)) involve summations whose index runs from 
1 to infinity. In practice, it is not possible to compute 
infinite sums. However, due to the fact that for all reason- 
able signals v(t), the magnitude of the coefficients gets 
smaller as k gets bigger, the infinite sum can be truncated 
without introducing significant computational error. Further- 
more, terms involving the magnitude squared of these coeffi- 
cients will tend to decay more rapidly allowing early truncation 
of the infinite sums. 

G. DISCUSSION 

In order to carry out the evaluation of Equation (2.90) 

to a reasonable degree of accuracy, the mathematical expression 

given by Equation (2.90) was evaluated by computer for fixed 

SNR values. Computation of P e was carried out for a fixed 

value on the upper limit on the index k. Recomputation of 

was carried out everytime this upper summation limit was 

incremented by 1. Recomputation of P^ was stopped after an 

incrementation of the upper limit in the sum terms did not 

yield an appreciably different value for P g in comparison to 

the P value before incrementation, 
e 

The result of these evaluations has been presented in 
Chapter IV in the form of tables and graphs while P g has been 
plotted as a function of SNR for different values of the 
parameter m. Recall that m specifies our relative prior 
knowledge about the signal's random time delay. 
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III. ALTERNATIVE SUBOPTIMUM RECEIVERS 



In the previous chapter, optimum detectors based on 
likelihood ratio tests were derived under conditions that are 
equivalent to low SNR conditions for signals having random 
time delay in the presence of additive white Gaussian noise. 
The tests derived gave rise to different receiver structures 
which, whenever possible, were analyzed in terms of their 
error probability performance. Both the receiver structures 
and their performance were seen to be a function of the uncer- 
tainty in the signal time delay or more precisely, a function 
of the p.d.f. of the signal time delay A. 

In this chapter, two suboptimum receivers are analyzed. 

The approach used here is different in the sense that re- 
ceivers are derived via strictly heuristic means and their 
performances are then analyzed. These receivers are: 

A. N-correlator Receiver 

B. Estimator-correlator Receiver 

A. N-CORRELATOR RECEIVER 

This receiver is shown in Fig. 3.1. The idea is to corre- 
late the received signal r(t) with N delayed reference signals 
generated at the receiver. If a signal is present, presumably 
at least one correlator branch would produce a large output so 
that the sum l could be guaranteed to be much larger under 
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Fig. 3.1. N-Correlator Receiver 

these conditions than when the correlators are fed by low 
level noise only. Thus, if l is greater than a prescribed 
threshold, then the presence of a signal with a random time 
delay is declared; otherwise only the presence of noise is 
declared. Here we assume for convenience that the observation 
time T q is an integral multiple of T, the period of the 
signal v(t). Observe furthermore that we are not required to 
assume low SNR conditions as done in Chapter II. 

In order to analyze this receiver in more detail, observe 

that 
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(3.1) 



i 



N 




T 

O n m 

/ r ( t ) v ( t 

0 



N 



l 

n=l 



r 

n 



where 



r 

n 



A 



T 

( ° nT 

/ r (t) v(t -^-)dt 



(3.2) 



If we consider the two hypotheses, as done in Chapter II, 
namely 



H x : r(t) = v(t-A) + n ( t ) 

H q : r(t) = n ( t ) 



0 < t < T (3.3) 

— — o 



then conditioned on hypothesis H Q , l is a Gaussian random 
variable (r.v.) with 




(3.4) 



since the noise is assumed to be zero mean, and 



E{r n r mlV = 



E{ / n(t)v(t-^dt / n(t) v(T-^-)dx} 



N o 

o ( nT. .. mT. 

— J v(t — N )v(t - lT )dt 



0 



N 

-£a 

2 nm 



(3.5) 
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where 



A 



nm 



A 



r ° nT, , , mT, ,, 

J v(t )v(t )dt 



are the signal discrete cross correlation terms. Thus, 



EU | H q } 



N 



I 

n=l 



E{r n |H 0 ) 



0 
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N N 

E{(T | H n } = E{ l l r r |H n } 

0 nil mil n 0 



N N N 

■ T I I 
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nm 



A 



c 



2 

i 



Now conditioned on both and X, £ is also a Gaussian 
r.v. with 

T o 

E{r | H. , X } = / v(t-X) v(t -^)dt = a (X) 

n 1 1 g J N n 

and 



( 3 . 6 ) 



( 3 . 7 ) 



( 3 . 8 ) 



( 3 . 9 ) 
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E{ [r n -E{r n |H 1 ,A}] [r m - E{ r m | H 1 , A } ] |H 1 ,A> 



(3.10) 
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so that 



N 

E{£|H ,A} = l a (A) 

n=l 
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Thus the conditional p.d.f.'s of £ are 



f L «l H o> 



2 2 

1 e " £ / 2CT £ 
V^TT of 
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a n ( A ) ] 2 } 
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(3.13) 
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N 



2 /_ 2 



f L ( £ | H L , A ) = 



- U - l a n (X ) ) /2a 



n=l 



-J 2tto 2 



(3.14) 



with 



f L ( £ | H 1 ) = / f (£ |H 1 ,A) f A U)dA 



(3.15) 



The detector of Fig. 3.1 performs the test 



£ 




H 



0 



(3.16) 



so that its error probability, P e is given by 



Y 



P e = P{H 0 } / P L U|H Q )d£ + P(H 1 > / p L (£|H 1 )d£ (3.17) 



Y 



Let p = P{H Q } and 1-p = P{H^} so that using equations 
(3.13) and (3.15), equation (3.17) becomes 



= P 



-- ~l 2 /2o 2 y - . 
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(3.18) 
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Interchanging the order of integration in the second 
(double) integral of equation (3.18), we obtain 
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(3.19) 



Thus, the receiver probability of error becomes 



P 

e 



N 



fy- l a (A) 
f u n 



p erfc* (~~) + (1-p) / f A ( A) erf *| 



n=l 



dA 



(3.20) 



The threshold y can be optimized to minimize . Such an 

approach unfortunately leads to a complicated integral equation 

that is not easily solvable. in practice however, P can be 

e 

calculated on the computer for y = 0 , and then incrementing 
| y | until an optimum value of y is found for a given SNR. 

It must be pointed out that the receiver of Fig. 3.1 will 
not always yield desirable results. For example, if v(t) is 
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sinusoidal, in the absence of noise, £ = 0 which means 
that with noise present, the receiver is unable to detect 
the presence of the signal. Fortunately for v(t) sinusoidal, 
a better receiver is available and its performance is well 
known [Ref . 4 ] . 

1 . Special Case of Triangular Wave 

Consider now a periodic waveform of triangular pulses 
and its two arbitrarily delayed versions as shown in Fig. 

3.2. If we represent one triangular pulse by p(t),then we 
can express v(t) as 



M 

v ( t ) = l p(t-jT) (3.21) 

j = 0 

In order to find the cross correlation terms A 

nm 

required for the computation of P g [equation (3.19)], 
instead of taking all the pulses in equation (3.21), we 
simply consider the cross correlation of a single pulse of 
period T with a delayed version of itself and then multiply 
the result by the total number of pulses M within the interval 
T , since v(t) is assumed to be periodic with MT = T q . 

In order to make computation simpler, Fig. 3.3 shows p(t) , 
p(t-a) and p(t-B), where a < B has been assumed. 
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v(t) 






Fig. 3.2 Periodic Triangular Wave and Two of Its 
Delayed Versions 
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Fig. 3.3 One Period Restricted Triangular Wave and 
Two of Its Delayed Versions 
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Observe that for a < 6 



/ v ( t-a ) v ( t-B ) dt = M 



a 



/ p[t+(T-a) ]p[t+(T-3) ]dt 



3 T 

+ / p ( t-a) p [t+ (T-B ) ] dt + / p(t-a)p(t-B)dt 

a 6 



= M 



a 



/ [At+A(T-a) ] [At+A(T-B) ] dt 



3 T 

+ / A (t-a) [At+A (T-B ) ] dt + / A (t-a) A ( t-B) dt 

a 3 



= ma 2 ( ^i + “|i - a 6 T + & - + ^ 



( 3 . 22 ) 



nT mT 

For a = — and 3 = equation (3.22) becomes 



ma 2 t 3 (JL- + JL-m+J^.JL+k) 

ivii\ i v 2 + 2N 2 + Z 2N + 3' 

2N N 2N 



= A 



nm 



(3.23) 



provided n < m. For n = m or equivalently a = 3/ the quantity 
inside the parenthesis in equation (3.23) becomes j. 

Similarly for a > 3, equation (3.22) becomes 

r ° 2 a^T aT^ 3^T BT^ 

/ v ( t-a) v ( t-B ) dt = MA (^ 2 ^-^i a g T + £_L +£i_ +__) (3.24) 
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(3.25) 



. n T mT 

Letting a = — and 6 = — in equation (3.24) , we obtain 



N 



A 



nm 



..,.2 3 . n' 

MA T ( 



2N 



n _ nm m m l. 

2N 2 + _ 2 + 2N 3’ 

N 2N 



provided n > m. Here again for n = m, the quantity inside 
the parenthesis in equation (3.24) becomes j. 

To proceed further with the analysis, the double 
summation of equation (3.8) is broken into three summations 
to account for n < m, n = m and n > m. Thus 



N 



1 

n=l 



N 



I 

m=l 



nm 



N N-l 

l l 

m=n+l n=l 



nm 



N 

+ l 

n=l 



nn 



N N-l 

+ I l 

n=m+l m=l 



A 



nm 



(3.26) 



Observe that the L.H.S. of equation (3.26) represents 
the sum of the elements of an N *N matrix. Whereas the middle 
term on the R.H.S. of this equation denotes the sum of the 
diagonal elements of the N xn matrix, the first term is the 
sum of all the elements in the upper triangular part and the 
third term is the sum of all the elements in the lower triangu- 
lar part of the N xN matrix. 

In order to evaluate the first term on the R.H.S. of 
equation (3.26), use of equation (3.23) (n < m case) yields 
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N-l N 

I 1 A 

n=l m=n+l 



= ma 2 t 3 l l { 



N-l N 



nm 



n=l m=n+l 2N 



n n _ nm m m l-, 

2 2N „ T 2 + „„2 2N 3^ 



N 2N 



2 3 
MA T J 



N-l N 



N-l N 



^ I. I . " 2 + 5TT I I " 



2n n=l m=n+l 



n=l m=n+l 



N-l N 



N-l N 



T l -l nm + — ^ l l 



T 

N n=l m=n+l 



m 



2N n=l m=n+l 



, N-l N , N-l N 

5 I ! »* I I 1 

n=l n=n+l n=l m=n+l 



HA 2 T 3 [ 3N - 4 2 »+ 1 ] 



(3.27) 



The second term on the R.H.S. of equation (3.26) 



becomes 



N 



1 A 



n=l 



nn 



MA 2 T 3 l i- 
n=l 3 



MA 2 T 3 [ j] 



(3.28) 



since A = — for n = m. 
ran 3 

Finally, substitution of equation (3.25) (n > m case) 
in the third term on the R.H.S. of equation (3.26) yields 



N N-l 

I I A 

n=m+l m=l 



nm 



MftV [-4 - W + K + W 4j 

2N 2 2N N 2 2N 2 2N 3 



MaV ; 3N -4N+1 , 



(3.29) 



62 



Adding the results of equations (3.27), (3.28) and 



(3.29) yields 



N 



1 

n=l 



N 



I 

m=l 



A 



nm 



M7 . 2 3 r 3N" 
MA T [ 



+ 1 
12 



] 



Thus from equations (3.8) and (3.30), 
variance of i becomes 



N N N 

~T l l 

n=l m=l 



A 



nm 






o „ ,3N +1, 

~t E ( — r ~ } 



since 



9 

j v (t)dt = M j v (t)dt 
0 0 



T 

M / (At) 2 dt 
0 

2 3 

MA T A „ 

T - " E 



is the energy of 
In order 
need to evaluate 
as 



v(t) for 0 < t < T . 

— — o 

to compute the probability of 
N 

J a (A) given by equations 
n=l 



N 



I 

n=l 



a n (A) 



N 



1 

n=l 



T 

r ° nT 

/ v(t-A)v(t — j^") dt 



(3.30) 



the conditional 



(3.31) 



(3.32) 



error P , we 
e 

(3.9) and (3.11) 



(3.33) 
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nT 

A Fourier series expansion of v(t-A) and v(t -— ) 



yields 



N 



n=l 



a n (A; 



. 2tt. , . £n T . 2tt 

3 .pkA - j 2 tt o ] (k+£)t 

L I I V k V £ e e / e dt 



N OO oo 



n=l k=-°° Z=-°° 



0 



• 2tt, , ... . ,2TTk. 

“ o -D-m-kA N Dn(— — ) 

T o . I l\l e l. e 

k=-« n=l 



(3.34; 



Observe that for k = ±iN, where i is an integer, the 
second summation involving the index n yields N, otherwise it 
yields zero. 

Thus, equation (3.33) can be written as 






cos ^-kNA) 



(3.35) 



since 



IV 2 - l v -kl 2 
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nT m.T 

In equation (3.6), if we represent v(t — — ) and v(t — — ) 
in terms of Fourier series, then 



N N 

y Y A 

L , L , nm 
n=l m=l 



N N °° °° 



I I I I W - e 

n=l m=l k=-°° £=-°° 



.2iTkn .27T£m T .2tt 

-3— -3— , ° 3-t 

J e dt 

0 



27Tk, 



. ,2iTk, 



= t y ivj"- y e ' 

o , L 1 k 1 L . 
k=-°° n=l 



2 v jn( — iT } ? ^ m( ~N _) 



l e 

m=l 



= T l | V, 

o , L 11 



k=-oo 



sin (kn ) 

• / k7r \ 
sin (*ir 



(3.36) 



since 



N 



n=l 



jn9 _ sin (N6/2 ) j(N+l)0/2 

J sin (0/2 ) e 



(3.37) 



Observe that the R.H.S. of equation (3.36) is always zero 
except for k = ±iN where i is an integer. Also observe that 
for k = ±iN, the R.H.S. of equation (3.36) becomes an indeter- 
minate form. Due to that reason, differentiating sin(kir) and 

k 7T 

sin(— ) w.r.t. k and then taking the limit k -* iN, equation 
(3.36) becomes 



N N 

l l A 

n=l m=l 



nm 



t n 2 y | V. .. | 2 

o , L 1 kN 1 

k=-oo 



t o n2(v o + 2 ziV 2 ) 

k=l 



(3.38) 



65 



so that equation (3.8) can be written as 



4 - + 2 1 W kN l 2 ) 

k=l 



(3.39) 



Equations (3.35) and (3.39) are of general form and can 
be applied to any signal, not just a triangular waveform. 

The Fourier coefficients of. the periodic triangular 

wave are 



V, 



. 1 



AT 

2 



( 3AT 
' 2irk 



k = 0 



otherwise 



(3.40) 



and 



„ 2 m 2 

A T 



k = 0 



V, 



2 2 
AT 

. 2 . 2 
4ir k 



otherwise 



(3.41) 



so that the use of equations (3.31) and (3.35) yields 



X _ I 

^ Z 3N^ + 1 

[N E (— — o~~) J 
o 8 



Y 



r E , 3N + 1, ! 1/2 

l N ' 8 

o 



(3.42) 



where y* = y/N 
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N 



I a (X ) 

L n n 
n=l 



8 



Iw-^t-t-7 

o 3N +1 



2 2 oo 

N(^~ + 2 l 



a 2 t 2 cos NX 



2 2 
k=l 4 tt N 



2 2 
AT 



+ 2 l 



2 2 
AT 



2 2 2 
k=l 4ir k N 



since 



8 



/— (■ 

N o 3N 2 +1 



NliH- 2 [l-i + (|) 2 ]} 

N 

1 + 1/m 2 



(3.43) 



I 

k=l 



cos k0 



"-iie+V o < e < 2 tt 

6 2 4 — — 



(3.44) 



and 



N 

Y " l ct 
n=l 



n 



(X) 



a 



£ 



/_£( 

r El ,3N 2 + 1 n1 1/2 V N o 3N 2 +1 

1 N 1 8 H 

o 



tj/i , 2rl X , ,Xv2-|'i 

U ? * " T + T 1} 
1 + 1/3N 2 



(3.45) 



Finally, the expression for P , after substituting equations 
(2.1), (3.42) and (3.44) in equation (3.20), becomes 
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= p erfc* 



f E ( 3N +1 ) j 1/2 



1+ (1-P) 



2ir A 

T/2 e m COS T ' 



N 



-T/2' 



TI q ( m) 



erf* 



r E (3N +1) i 1/2 

1 N 8 J 

o 



/-£(_§. 

(N ' ,..2, , 
o 3N +1 



-) 



3N 3 +6N [i + (^) 2 ] 



(3N 2 +1) 



dA 



(3.46) 



As mentioned earlier, the receiver of Fig. 3.1 has 

been designed via ad-hoc methods, so that the optimum threshold 

setting is not readily obtainable. The difficulty in applying 

an analytical approach to the problem of finding the optimum 

dP 



threshold arises in the solution of 



d Y 



= 0 for y • Since 



N 

" (Y- l a n (A) ) 
n=l 



dP 



d Y 



e 1 

= "P — e 

/2 7T 0 „ 



- %2/2a \ , T/2 1 

1 + (1-p) / f.(A) —i- 



2a! 



-T/2 



/2ttu , 



dA 

(3.47) 



then, solution of 



dP 

e 

d Y 



= 0 is equivalent to solving 



N 



N 



-T/2 



T/2 

/ 



Y l a ( A) t I a r, (^) ] 
n=l n n=l n 

2 o 2 

a „ 2a„ 



(A) e 



dA = 



(1-p) 



(3.48) 



68 



Clearly, the computer could be used to search for a value of 
Y that satisfies the equality prescribed in equation (3.48). 

The approach that has been used in trying to obtain 
the optimum threshold involves using the computer to evaluate 
equation (3.46) as a function of y for fixed SNR. Thus the 
value of y for which P g is minimum yields the optimum thres- 
hold for that particular value of SNR. This procedure is 
repeated for different values of SNR. The resulting values 
are then used to plot P e vs. SNR. These curves have been 
plotted using optimum threshold settings and the results have 
been summarized and discussed in Chapter IV. (Observe that 
the threshold y is SNR dependent. This situation is typical 
of noncoherent signal detection problems.) 

B. ESTIMATOR-CORRELATOR RECEIVER 

The second suboptimum receiver being proposed is the re- 
ceiver structure shown in Fig. 3.4. it is basically an 
estimator-correlator receiver in which first a coarse estimate 
of the time delay in v(t) is made and then this information 
is used to process the delayed signal r(t) coherently with 
an estimated reference. 

In this receiver the largest output of the N correlator 
branches is used as an estimate of the signal phase by the 
logic choice. This estimate is used to provide an accurate 
local reference with which to perform the correlation detection 
operation. The output £ of the main correlator is used for 
threshold comparison in order to make the decisions. 



69 



Delay 
To sec . 



4 



f 0 
U1 U ul 
QJ 0-h 







70 . 



Fig. 3.4 Noncoherent Correlator Receiver 



Thus 



T 

O 

£ = j r(t)v(t-a)dt (3.49) 

0 



T 2T 3T 

where a = . . . , T and a may take on only one such 

value for 0 £ t <_ T q . In order to determine receiver perfor- 
mance, we note that £ is a conditionally Gaussian random 



variable with 



E{ £ | H Q } = 0 (3.50) 

T 

o - 

Var{£|H Q } = E{ [ / n ( t) v ( t-a) dt] } 




2 



v 



(t-a) dt 



N OO oo 

= ~Y l l 

k = -°° £ = -00 



V k V £ 



“j^r(k+£)a 



T . 2 tt ,, , „ . . 
o : — (k+£)t 

/ e dt 

0 



N T 
o o 

2 





(3.51) 



where are the exponential Fourier series coefficients of 
v(t) , and E is the energy of v(t) for 0 £ t £ T. 
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Also 



T 

o 

E{£|H,,A} = / v ( t-A ) v ( t-a ) dt 

0 



l l 

k=-oo £=-00 



• -j^(kA+£a) 

W 




j±f(k+£)t 
e dt 




-j— k(A-a) 
e 



= T o [ I Vq I 2 + 2 l |v k | 2 cos^k(A-a)] = m(A) 

k — 1 



(3.52) 



It can be shown without a great deal of difficulty that 



Var{ £ | , A } = o 2 



(3.53) 



Therefore the conditional p.d.f.'s of £ are given by 



fiUlv 



i e - z2/2o l 

/27a t 



(3.54) 



, - [£-m (A) ] 2 /2c 2 

f T ( £ | H. , A ) = — - — e 

" 1 ^£ 



(3.55) 
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so that 



- -[£-m(A)] 2 /2 a 2 

f T ( £ | H. ) = / — - — e x f . ( A ) dA 

L 1 -co' / 27 o £ A 

Thus the expression for is given by 
P e = P{£ >y|H 0 )P{H 0 } + P{ £ <y|H 1 )P{H 1 } 

0 ° , -£ 2 / 2 a 2 Y 00 1 

= p / e Ji d£ + (1-p) / / _ c 

Y / 2 tto^ — oo — .oo ✓ 2*0 

X f A ( A ) dAd£ 

where 



p = P{H Q } and (1-p) = P(H 1 ) 



as before. 

After a change of variables, equation (3.57) can 
written as 



P 

e 



°° . 2 oo 

P / e X ^ + (1-p) / 

Y/ 0 ^ /2tt — OO 



/ 



(Y-m(A) ]/o Q 



/2lT 



= p erf c* (- 7 -) + ( 1 -p) / 



f A ( A) 



erf*( Y ~ m(A) )dA 

°£ 



(3.56) 



-m ( A ) ) 2 /2 o 2 



(3.57) 



be 



x / 2 dxf A (A)dA 



(3.58) 
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Since the receiver was not designed from likelihood ratio 
test principles, the threshold y must be set for minimum p . 
The approach used to find the optimum threshold involves 



dP 



0 

solving = 0 for y as before. 

Since 



J_ _ 



y 



/n e/2 

o 



y’ 



/E/2N 



where y* = y/N 



(3.59) 



and 



m(A) 



= /2E/N 



v 0 I 2 + 2 k ^ ] _^ V k ^ 2 

|V 0 I 2 + 2 l |V I 2 
U 1=1 ^ 



(3.60) 



we have 



y-m(A) _ 



T - /2E/N 
/E/2N 



v 0 | 2+2 I |V | 2 cos^(X-a) 
U k=l K 

00 

I V I 2 + 2 l I V | 2 
U 1=1 36 



( 3 . 61 ) 



Thus equation (3.56) becomes 



2ttA 

T/2 mCOS T 

p erf c* ( Y - ) + (1-p) / 6 



/ E/2N 



-T/2 



TI q ( m) 



x erf, 



Y /2E/N_ 

/E/2N ° 



V ° ' 2+2 J-l |V k' 2cOS ^ i(X " a) 

! v 0 ! 2 + 2 I l v £ l 2 
0 £=1 36 



dA 



(3.62) 
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As before, determination of the optimum threshold by solution 
dP 

0 

of ^ = 0 involves solving an equation of a form similar to 

that of equation (3.47). 

In general, solutions in closed form do not appear tracta- 
ble. A search for an optimum y may be performed or solutions 
for specific cases may be possible. This is illustrated by 
considering two special cases in which v(_t) is a 

1. sine wave 

2. square wave. 

1 . Sine Wave 

As a first special case, let 



v(t) 



A cos 





so that 



V 



k 




|k| = 1 

otherwise 



(3.63) 



(3.64) 



and 



oo 



2 l 

k=l 





(3.65) 
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Thus equation (3.60) becomes 



m ( A ) 



= /2E/N 



cos ^-(X-a) 

A 2 /2 



2tt 



/2E/N q cos ^-(A-a) 



(3. 6.6) 



and for p = (1-p) = j, equation (3.58) becomes 



. il 



erfc* (- 



Y 1 



T/2 2itA 

' mcos— — 



/E/2N 



-) + 



TI 0 (m) -T/2 



/ e 



x erf*(- 



Y : 



2tt 



/E/2N 



/2E/N q cos-^-(A-a) )dA' 



( 3 . 67 ) 



As mentioned earlier in Chapter II, m = 0 corresponds 
to a uniform p.d.f. on A, i.e. p(A) = i, so that the error 
probability expression for m = 0 becomes 



' erfc* ( 

( 



Y 



/E/2N 



T/2 

: ) + ;jr / erf * (- 
i -T/2 



r- 



2 7T 



/E/2N 



- /2E/N cos^-(A-a) ) dA 



(3.6 8) 



2 . Square Wave 

Consider now the second special case in which v(t) 
is a periodic square wave as shown in Fig. 2.6. In Chapter II 
we showed that 

for k = odd 
for k = even 
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so that equation (3.59) becomes 



m(A) 



= /2E/N 



00 cos k2ir(^7fr^) 

1 2 “ 

k=l k z 



00 

i 4 

£=1 £ 



= /2E/N 



tt 2 (| - (A^L) + (^) 2 ) 

' u 2 /6 



= /2E/N o (l - 6(^L) + 6(^) 2 ) 



(3.69) 



Therefore 



P 

e 



1 

2 



erfc* (- 



Y 



/E/2N 



+ 



1 

(m)T 




-T/2 



mcos 

e 



2nX 

T 



erf* ( - — - /2E/N (1 - 6(^=2-) + 6(4^) 2 )d£( 

/E/2N q o T T ) 



(3. 70) 



Observe that equations (3.35) and (3.6 2) involve 
summations whose indices run from 1 to infinity. In practice, 
it is not possible to compute infinite sums. However, due to 
the fact that for all reasonable signals v(t), the magnitude 
of the coefficients gets smaller as k gets bigger, this 
infinite sum can be truncated without introducing significant 
computational error. 

Fortunately for triangular and square wave signals 
considered as special cases in this chapter, the argument of 
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summation in equations (3.35) and (3.61) turned out to be of 
such a form that it could be expressed in closed form. 

For the sine wave, summation disappeared due to the fact 
that its coefficients exist only for k = ±1. These results 
in the form of graphs are presented and discussed in Chapter 
IV. 
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IV. DISCUSSION OF GRAPHICAL RESULTS 



This chapter presents graphical results obtained by 
applying the analytical results of the previous chapters to 
specific examples and carrying out the required computations 
on the computer. The plots are intended to display receiver 
performance in terms of probability of error (P ) as a function 
of SNR. Some curves of P g versus comparator threshold level 
have also been included in order to show the dependency of P e 
on the threshold which in turn depends on SNR. The graphical 
results displaying performance of the optimum receivers are 
presented first, followed by graphical results displaying 
performance of the suboptimum receivers analyzed. 

A. GRAPHICAL RESULTS FOR OPTIMUM RECEIVERS 

1 . Receivers Operating on Signals with Nonzero DC 
Component (Uniform p.d.f. on A) 

In Chapter II, using statistical communication theoretic 
principles, a simple integrator receiver (Fig. 2.2) was derived 
that could discriminate between signal plus noise and noise 
only hypotheses provided the signal had a non-zero DC component. 
Recall that the receiver derivation and the subsequent per- 
formance evaluation carried out was made possible by the 
two-term approximation on the exponential appearing in Equation 
(2.19). Equivalently, low SNR conditions were assumed so that 

the P results are only valid for low SNR's, 
e 
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Figure 4.1 Performance (Uniform p.d.f.) 
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Figure 4.1 shows the performance of this receiver as 
a function of SNR. This graph was obtained through numerical 
evaluation of Equation (2.40) under the assumption that v(t) 
has a strong DC component in comparison to its harmonics, 
equal prior probabilities, and uniform p.d.f. on the r.v. A. Due 
to the constraint imposed by Equation (.2.20) , Equation (2.40) 
is valid up to an SNR of about -4.1 dB . Thus Fig. 4.1 is 
plotted only up to the SNR of -4.1 dB . . 

Observe that P g is at a very high level of about 0.3. 
However, it decreases as SNR increases, as expected. While 
it is clear that no receiver could operate with an error 
probability of 0.3, the figure does give an indication of the 
performance level that can be expected at low SNR's. For 
high SNR's, it is expected that P g will continue to decrease 
to reasonable levels. However, for high SNR's, the receiver 
of Fig. 2.2 is no longer optimum. 

2 . Receiver Operating on Arbitrary Signals (Nonuniform 
p.d.f. on A ) 

The two-term approximation on the exponential appearing 
in Equation (2.19) and application of the nonuniform p.d.f. on 
A given by Equation (2.1) resulted in the receiver of Fig. 2.4 
that could operate on signals having either zero or non-zero 
DC component. Observe that as m -»■ 0 , the receiver of Fig. 

2.4 becomes equal to the receiver of Fig. 2.2. Two specific 
signals, a sine wave and square wave, both having zero DC 
component, were analyzed. 



81 



Figure 4.2, obtained through numerical evaluation 
of Equation (2.80), shows the performance of this receiver 
as a function of SNR when v(t) is a sine wave. Three differ- 
ent values of the independent parameter m were considered. 
Recall that m controls the level of uncertainty about the time 
delay A. As expected, as m increases, P^ decreases for a 
given SNR. For example, ai SNR of -5.2 dB is required in 

order to achieve aP of 0.4 at in = 2. With m = 90, we need 

e 

an SNR of -9 dB to achieve the same P e . Note that the m = 0 
case has not been considered here due to the fact that m = 0 
corresponds to a uniform p.d.f. on A. The receiver for that 
case (Fig. 2.4) is unable to detect signals having zero DC 
component, resulting in P g equal to 0.5. 

Curves shown in Fig. 4.3 are the result of numerical 
evaluation of Equation (2.90) in which three different values 
of the independent parameter m have been considered. The 
m = 0 case has not been considered here either due to the 
same reasons presented for the sine wave case. 

Observe that Equation (2.90) involves the computation 
of infinite sums which can not be carried out in practice. 
However, due to the fact that for all finite power signals 
v(t) , the magnitude of the coefficients gets smaller as k 
gets bigger and thus, truncation is possible. Inclusion of 
more coefficients in Equation (2.90) corresponds to increased 
accuracy in the results. Furthermore, P g changes as k is 
incremented, however up to a point only as further increase in 
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Figure 4.2 Performance (Sine Wave) 
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Figure 4.3 Performance (Square Wave} 
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k beyond this limit do not result in an appreciable change 
in P . This phenomenon is illustrated in Table I. 





table I 




probability of 


ERROR VS tt OF 


COFFICIENTS SUMMED 


( 


SNR = -10 OB) 




* OF COFFICIENTS 


PROBABILITY of 


SUMMED 




ERROR 


I 




0.4451 24370+00 


2 




0. 4 A A3 3 947 0+ 00 


3 




0 .444325500 + 00 


4 




0.444325430+00 


5 




0.444325420+00 


6 




0.444325420+00 


7 




0.444325420+00 


a 




0.444325420+00 


9 




0.444325420+00 


1 0 




0.444325420+00 



We see that keeps on decreasing until the number of 
coefficients summed is 5. For k > 5, P g remains unchanged. 
Thus, for this example the computation of Equation (2.90) was 
carried out with sums truncated to 5 terms. The performance 
curves for the square wave signal are shown in Fig. 4.3 
for low SNR conditions. With m = 2, aa SNR of -4.1 dB is 
required in order to achieve a P g of 0.4 while for the same 
P e , with m = 90, the SNR required is -8.8 dB . 

B. GRAPHICAL RESULTS FOR SUBOPTIMUM RECEIVERS 

In Chapter III, two suboptimum receivers were analyzed. 

The design approach used in that chapter was different in the 
sense that receivers were derived via strictly heuristic 
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means. Therefore, the optimum threshold setting for the 
receiver was not readily obtainable. However, a computational 
approach was used to solve the problem of properly setting 
the threshold prior to evaluating P g as a function of SNR. 

1 . N-Correlator Receiver 

The first suboptimum receiver analyzed in Chapter III 
was the N-Correlator receiver of Fig. 3.1, and its performance 
with v(t) a periodic triangular pulse was evaluated. The 
graphical results corresponding to the numerical evaluation 
of Equation (3.46) are shown in Figs. 4.4 through 4.12. As 
mentioned earlier, using the computer. Equation (3.46) has 
been evaluated as a function of the threshold y for fixed SNR 
values starting at -10 dB and increasing in increments of 5 
dB . For a given SNR, the threshold value for which P is 
minimum corresponds to the optimum threshold for that particu- 
lar SNR. . These minimum points corresponding to each SNR have 
been used as optimum threshold values in Equation (3.46) and 
the resulting P g has been plotted in Fig. 4.12 as P g vs SNR. 

Observe that Equation (3.46) involves the previously 
encountered parameter m and N, namely the number of correlators 
used by the receiver. Two different values of each parameter 
m = 0, 90 and N = 2, 16 have been used in the computations 
in order to evaluate their effect on receiver performance. An 
inspection of Fig. 4.4 (m = 0 , N = 16) and Fig. 4.6 (m = 30, 

N = 16) reveal that for a given N, m has no significant effect 
on p . On the other hand Fig. 4.4 (m = 0 , N = 16) and Fig. 

4.8 (m = 0 , N = 2) reveal a small effect of N on P g for a 
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given value of m. The actual effect of N on P is shown in 

e 

Fig. 4.12 where has been plotted as a function of SNR for 
two values of N. Observe that the higher value of N yields 
better performance. To achieve a P^ of 10 -4 with N = 2, the 
SNR required is ~ 13.5 dB whereas for the same P g with N = 16 
we need an SNR ~ 13.2 dB . This result is once again expected 
as using few correlators might result in a small output on 
the correlator branches when v(t) received in r(t) has a 
delay very much different than the locally generated v(t) . 

While it is possible to "see" the reason for the 
independence on m in Equation (3.46), a heuristic explanation 
of this observed result follows the following argument. The 
receiver of Fig. 3.1 correlates the incoming signals with local 
replicas of v(t) delayed by multiples of T/N. Depending on 
the actual delay in v(t), upon reception of r(t), one of the 
correlator branches will produce the largest output (or there 
may be a tie between two branches) . The output produced, 
however, namely which correlator branch is largest, should 
be independent of the level of uncertainty about A. Clearly, 
as m increases, it is easier to predict which correlator 
branch output will be largest. But this is inconsequential 
as the receiver is only trying to detect the presence of v(t) . 

2 . Estimator-Correlator Receiver 

The second suboptimum receiver analyzed in Chapter III 
was the Estimator-Correlator of Fig. 3.4. Its performance 
for v(t) either a periodic sine wave or a square wave was 
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Figure 4.5 P e vs Thld (m 
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Figure 4 . 8 
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Figure 4.9 
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Figure 4.10 P e vs Thld (m = 30, N = 2) , (SNR = 
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Figure 4.12 



considered as a function of SNR. As in the case for the N- 

correlator receiver, this receiver does not have an optimum 

threshold that is obtained from the theoretical development. 

However, a procedure similar to that used for the N-correlator 

receiver has been used to obtain the optimum threshold setting. 

The P has been evaluated as a function of SNR except that 
e c 

here, plots showing P g as a function of y have not been included. 

Figures 4.13 through 4.16 have been obtained through 
numerical evaluation of Equations (3.67) and (3.70). Two 
curves in each figure have been plotted to show the effect 
of the parameter m on the receiver performance. Observe that 
Equations (3.67) and (3.69) involve a parameter a which can 
only take on values iT/N, where i = 1,2,...,N. Observe that 
i = N corresponds to a = T or equivalently a = 0 since v(t) 
is T periodic. 

Two values of a have been considered for each of the 
signals, namely the sine wave and square wave, in order to 
determine its effect on receiver performance. An inspection 
of Fig. 4.13 and Fig. 4.14 displaying the performance of the 
receiver for v(t) a sine wave, or inspection of Fig. 4.15 and 
Fig. 4.16 displaying the performance of the receiver for v(t) 
a square wave, it can be observed that a does not have any 
significant effect on P g for either case. The logic support- 
ing this observation on receiver performance as a function of 
a follows from the fact that a is just the estimate of the 
time delay of the signal, which is obtained by choosing that 
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correlator branch that produces the largest output. The 
actual value of a is irrelevant in so far as the performance 
of the receiver is concerned. The receiver simply uses the 
estimate of the time delay to produce a decision based on a 
correlation operation. However, one of the factors that might 
improve or degrade the error probability is the accuracy of 
the estimate of the time delay of the signal, which is 
itself a function of N and SNR. 
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Figure 4.13 Performance (Since Wave)., 
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Figure 4.15 performance (.Square Wave) , 
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Figure 4.16 Performance (Square Wave), (Alpha 



T/21 



102 



V. CONCLUSIONS 



In this thesis, the analysis and design of noncoherent 
digital receivers that optimally (and in some cases subopti- 
mally) detect an arbitrarily shaped periodic waveform with 
random time delay, has been carried out. The probability 
density function (p .d . f . ) of the random time delay A is assumed 
to be of some known form. A uniform p.d.f. and a (parameter 
varying) non-uniform p.d.f. have been considered and their 
effect on the receiver performance has been studied. 

The optimum (in minimum probability of error, P , sense) 
receivers have been designed on the basis of likelihood ratio 
test under low SNR assumptions. The design of suboptimum 
receivers is based on heuristic approaches that intuitively 
yield structures or algorithms that yield good performance. 

The analysis carried out in this thesis, allowed for the 
possibility of unequal prior probabilities of occurrence of 
the hypotheses, however all examples and computations assumed 
equal prior probabilities for simplicity. 

The low SNR assumption allowed certain mathematical 
complexities to be overcome. Therefore graphical results 
involving the performance (P ) of optimum receivers are valid 
and therefore presented only for SNR values less than or equal 
to -4 dB . Operating at SNR values less than or eciual to 0 dB 
is the basic requirement of covert communication and may 
also be necessary in certain radar applications. 
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The optimum receiver operating on the signals that have 
non-zero DC component, detects the signal with higher level 
of error probability and it decreases with the increase of 
SNR as expected. While it is clear that no receiver could 
operate with such an undesirable level of error probability, 
the graph displaying the performance (P ) does give an indi- 
cation of the performance level that can be expected at low 
SNR's. In the region of high SNR, it is expected that error 
probability will continue to decrease to reasonable levels. 
However, for high SNR's this receiver is no longer optimum. 

The second optimum receiver operating on arbitrary sig- 
nals does also yield poor performance. But here again we 
must not forget that the receiver is operating in the low 
SNR region. The graphs displaying the performance of this 
receiver for the sine wave and square wave give an indication 
of the performance that can be expected at low SNR's. Also, 
they reveal an improvement in the performance with an in- 
crease in the parameter m. It is concluded that for high 
SNR's, this receiver is also no longer optimum. 

One of the receivers designed via heuristic means was a 
(suboptimum) N-correlator receiver. Its performance was 
derived and tested for the case in which v(t) was a periodic 
triangular waveform. While this receiver is unable to detect 
a sine wave, it is capable of detecting a square wave or some 
other periodic waveform. The other receiver designed via 
heuristic means, namely the Estimator-Correlator receiver, 
was analyzed in terms of its P g as a function of SNR. 
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From the graphical results of the performance of the N- 
correlator receiver, it is concluded that receiver performance 
gets better as N, the number of correlators, gets bigger. 

It is further concluded that m, the parameter that determines 
the uncertainty about the time delay A, does not have any 
effect on the receiver performance. In general, it was 
observed that except for the case of the N-correlator receiver, 
the parameter m had a significant effect on the performance 
of the receivers. 

While a number of practical problems have been addressed 
in this thesis, related problems and potential studies abound. 
For example, the receiver obtained using a three-term approxi- 
mation in the exponential of Equation (2.19) needs to be 
analyzed further in terms of performance for both uniform and 
non-uniform p.d.f. on A. Furthermore, other (suboptimum) 
receivers need to be proposed and analyzed, such as simple 
energy detection receivers. 
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APPENDIX A 



DERIVATION OF MEAN AND VARIANCE OF R.V. A 



The p.d.f. of r.v. A (time delay) given by Equation (2.1) 



is 



2 TT 

m cos— A 

f A ( A ; m) = T t ( m ) 



A < I T/2 I . 0 < m < 



(A.i; 



Substitution of Equation (2.41) into Equation (A.I) yields 



f A ( A ;m) 



= P 



l I D (m) COSp^A 

= — oo ^ 



TI Q (m) 



A < I T/2 I . 0 < m < oo (A. 2] 



so that the expected value or mean of r.v. A can be expressed 
as 



T/2 

E{ A } = / A f . ( A ; m ) d A 

-T/2 yv 



TI o (m) P =- 



J I_(m) 



T/2 

/ cos (p2r A/T) dA 

-T/2 



= 0 



(A. 3) 
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since the product of X (odd function) and cos(p2TfX/T) (even 
function) is an odd function and integral of an odd function 
produces zero. 

The variance of X is given by 



E { X } 2 



T/2 



-T/2 



X 2 f A ( X ;m) 



dX 



1 

TI 0 (m) 



T/2 



-T/2 



.2 mcos27rX/T ,, 
X e ' dX 



(A. 3) 



since we have a zero mean random variable. 

The variance is clearly a function of m. However, the 
integral of Equation (A. 3) is not readily evaluated in closed 
form. The computer has been used to evaluate 



Var { X } 



4l rt (m) 



0 ' ' 0 



r 2 mcosTTx , 

J x e dx 



(A. 4) 



and it has been clearly demonstrated that Var{X}/T becomes 
smaller as m grows. 
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APPENDIX B 



INTERPRETATION OF A SMALL EXPONENTIAL 
IN EQUATION (2.22) AND EQUATION (2.60) 

In Equation (2.22) and Equation (2.60) it has been assumed 

that 

T 

2 0 A 

rj— / r(t)v(t-A)dt = £ << 1 (B.l) 

o 0 

This appendix is devoted to analyzing under what conditions is 
this assumption valid. Observe that e is a r.v. so that 
Equation (B.l) must be analyzed in a statistical sense. Recall 
that 



E{ e | H Q } = 0 (B . 2 ) 

since the noise has been assumed to be zero mean and 

T 

EUlf^} = / v 2 ( t- A ) = j^-E (B.3) 

o 0 o 

where E is the signal energy for 0 £ t £ T^ . 

Thus, provided E/N q << 1, e under either hypothesis on 
the average takes the value zero or a value close to it. The 
excursions of e away from the mean can be obtained by evaluat- 
ing the conditional variances. Since 
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E{ e 2 | H 



and 



E{e 2 | H 



Thus 



- 



n ( t) v ( t-A ) dt] 



T 

4 ,, 0 

= // E{n ( t ) n (t ) } v ( t-A ) v ( t-A ) dt dx 

N 0 
o 



= t -I v (t - x)at = ST e 

O 0 O 



L } = E{t ~ / v 



T T 

o 0 p o 0 

M - / v ( t- A ) dt + 7 ^-J n ( t) v (t-A ) dt] 

o 0 o 0 



?F 9 4 , u 

E{(^) Z + 2 (-4 / n ( t) v (t-A) dt) 



2 

N 0 
o 



+ (— / n (t) v (t-A) ) 2 } 

o 0 



,2E, 2 2E 

N 

o o 



Var{ e | H Q } 



ST = Var{e l H i } 

o 



(B . 4 ) 



} 



(B . 5 ) 



(B . 6 ) 
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Therefore, provided that E/N 0 << 1, e remains very close 
to its average value which is either zero or very close to 
zero. Note also that in a more rigorous sense, we can bound 



P{e 2 >6 | H ± > < 




(B . 7 ) 



If we set 6=1, then 



p{ £ 2 >i|h 0 } < 



2E 

N 

o 



and 



(B . 8 ) 



PU 2 >1|H} < (|^) 2 + P (B . 9 ) 

o o 

In both cases if 2E/N << 1, then p{e > 1 or e <-l} is also 

o — — 

small. Therefore, a small E/N Q guarantees that the assumption 
of a small exponential in Equations (2.22) and (2.60) is 
valid with a high degree of probability. 
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